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In this paper we examine the possibility, considered 
previously by others, that the recently observed pulsed radio 
sources are pulsating white dwarfs. We first present the results of 
radial mode eigenvalue calculations (for various central densities 
and compositions) which were obtained by integrating the eigenvalue 
equations themselves rather than by using variational techniques. 

Computation of non-linear white dwarf pulsations are then 
presented, in which a hydrodynamic code was used to follow 
finite amplitude pulsations driven by pycnonuclear reactions. 

Very large amplitude pulsation of neutron stars are also 

ft 

considered. 

The possibility that the optical pulsation periods are 

« 

twice as long as the radio periods has arisen from observations. * 


We discuss a method for such a "frecruencv doublinq" based partly 
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The following white dwarf eigenf requency calculations were 
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obtained by computer integration of the Newtonian equations for 
the eigenfunctions. The eigenvalues were obtained by finding 
those values of the period for which the eigenfunction 5 = 6r/r 
and d§/dr are finite at the outer boundary and match to an analytic 
expansion (about the surface) of the eigenfunction. The eigen- 
values compare very well with those determined by variational 
2 

methods , but the method described here has the advantage of 

giving accurate eigenfunctions without having to postulate 

polynomial forms. In table 1 we have used the degenerate electron 

equation of state as first applied to white dwarfs by Chandrasekhar^ 

The effect of Hamada-Salpeter corrections to the equation of 
4 

state , etc., is found to reduce the periods slightly. In table 2 

we show the results of calculations corrected for the Coulomb 

Thomas-Fermi and exchange interactions (Hamada-Salpeter corrections) 

In both tables the ratio of nucleons to electrons was taken to be 

2.0. Thus the un -corrected models can represent for example 
1. 2 X 2 4 

C x or C + He compositions while the corrected model assumes 

12 

Z = 6 as with C 

For higher densities than those considered here, we have 

5 

found that relativistic effects are too large to neglect , 
and that periods considerably lower than 1.64 sec are very 
difficult to obtain. This difficulty, however, can be circumvented. 
We later give a mechanism which generates two radio frequency 
pulses for each stellar pulsation. Also, in the following sections 



we discuss dynamical calculations with simulated nuclear excitations, 
which tend to show that pulsation periods shorter than the 
fundamental may in any case be obtained. Thus there are at least 
two possible ways of producing models with radio periods shorter 
than the fundamental periods. 


Dynamical Calculations 

The full, non-linear hydro-dynamical equations of motion 
were used to follow the oscillations of a white dwarf starting 
from initial conditions close to equilibrium. We employ the 
Lagrangian form^ 

^7 = -4nr^ m (P+Q) - 


where m is the mass interior to radius r = r(m,t) , p is the 
pressure, and Q is a viscosity term introduced both to spread 
out shock fronts over several mass zones, and to eliminate 
the effects of the initial conditions on the final state of motion 

of the white dwarf. We employed 20 mass zones of equal mass, and 

7 , 8 

used a form for Q given by Colgate and White and Arnett as 
Q k+H = Cp k+is (v k+ls -v k ) 2 for v k+1 < v k and = 0 etherise, where 

v^ is the velocity of k-th mass shell, ^h e density, and 

Q 

C is a constant. Arnett chooses C «= 1 , but we chose C = 10 in 
order to damp out the effects of the initial conditions in 


reasonable computer time. 



The pressure is given by two terms P = P + p , where P 

e n ^ 

is the usual degenerate electron pressure and P is due to 

n 

9 2 

pycnonuclear reactions . We have used the form P n (p - p 1 ) 

where is a density at which the excess pressure from the heating 

from the nuclear reactions was assumed to die away. The reactions 

were assumed to cut on at a density p 2 > P^- The total mass in 

32 

the model used was assumed to be M = B.38 x 10 g, corresponding 

6 3 

to a central density P c of about 2 x 10 g/cm . Of course, no 
pycnonuclear reactions actually occur at this density. We would 
have liked to have run models with densities high enough for 
these reactions to occur, but we found that for P c ^ 10 gm/cm' 3 , 
the initial deviations from equilibrium of our model gave rise 
to gravitational collapse! However, one can imagine that the 
general results obtained would be applicable to more compact 
configurations which will not collapse for sufficiently small 
oscillations. We conclude, however, that large amplitude 
oscillations would result in collapse for the more compact 
configurations, at least down to the neutron star stage, as 
discussed below. 

We have assumed, that the nuclear energy source is very close 

4 12 

to the surface, where for example He might be converting to C 

2 

Using the form P = C (p - p _ ) for the over-pressure in the next- 

n o 1 

11 4 

to-last mass zone, we found that C Q = 9.0 x 10 and p^ = 5.0 x 10 
give reasonable oscillations with P 2 = 7.0 x 10- 1 - 4 . Figure 1 
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s’nows the final asymptotic time dependence of the outside radius 

and the central density. The time intervals between the large 

radius . maxima are nearly equal to the fundamental frequency for 

small pulsations. Note the appearance, however, of the small 

secondary maxima. In more realistic models, the secondary maxima 

might well be larger, thus shortening the apparent 'oscillation 

period. Many optical variable stars have 1 been noted to have 

periods shorter than the fundamental^ 

Similar calculations were also carried out for very large 

amplitude degenerate neutron star pulsations. For a neutron 

33 

star of total mass M ■ 1.0 x 10 gm, an oscillation period of 

1.33 sec can be obtained if the maximum radius is taken to be 

« 

2600 km. In Figure 2 we have graphed the outside radius as a 
function of time. Note that the period agrees with the estimates 
of Hoyle and Narlikar^. However, even if their mechanism to 
avoid high neutrino energy loss at maximum contraction is possible, 
we find that the density at ..iaximum extension is so low that free 
neutron decay takes place again resulting in neutrino loss. In 
fact, the damping time for the pulsations is of the order of a 
few months. 

/ 

Radio Emission 


In this section we discuss plasma mechanisms for generating 
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the observed radio emissions. We assume that the basic energizing 

mechanism is an electric field in a layer where the scale height 

is smaller than the mean free path of thermal electrons. Thus 

we have a "low density" plasma in which ordinary MHD is invalid, 

and the electric field may have a component parallel to the 
12 

magnetic field . From Faraday's law, -B = c7xE , we see chat 

maximum E is obtained for maximum | b| , which corresponds to maximum 

| dp/dt| in the deeper layers of the star, where MHD is valid. We 

will see that fields B > <>0 gauss are required, and thus the 

magnetic pressure dominates the plasma pressure in the low density 

layer. Changes in B in the interior will then be communicated to 

the atmosphere. But max jdp/dtj occurs twice during a pulsation 

cycle, and thus maximum E occurs twice during a cycle, making the 

stellar pulsation period twice the radio period. 

The appearance of an electric field E with a component 

parallel to the magnetic field will then give rise to longitudinal 

plasma oscillations in which ■‘-he electrons move parallel to B, 

at the plasma frequency u> . These oscillations themselves do not 

radiate energy since there is no magnetic field associated with 

them, but in an inhomogeneous medium there is coupling with 

13 

other radiative plasma wave modes , and the emitted radiation 

14 

would then be at the plasma frequency w . Kundu and Chitre 

P 

have also mentioned the possibility of plasma oscillations as 


an exciting mechanism. 
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Let us discuss radiation at 100 Me. The corresponding plasma 

8 — 1 

frequency is 2rr x 10 sec , and the electron density is then 

n = 1.24 x 10 8 cm . Let us assume that the thermal temperature 

is 10^ °K, and that the local density scale height is 1 0 = 6 x 
3 

10 cm, a typical scale height for a white dwarf atmosphere. The 

15 

electron mean free path \ is then , with In A = 15 (where In 
A is the Coulomb logarithm) , 

5 2, 3 

X e = 1 . 3 x 10 T /n In A = 7 . 2 x 10 cm. 

Thus, \ e > 1 0 , and MHD breaks down. The 1 0 assumed here may 
actually be too small, since the intense radio output from the 
pulsar will probably distend the atmosphere. This same effect, 
however, will also raise the temperature, and we will still have 
X e > 1 0 . 

We assume that the average power emitted from this layer is 

2p 16 

about 10 erg/sec. The conversion from longitudinal 

oscillations to radiative plasma waves could occur in many *(say 
7 

10 ) bursts per fundamental pulsation, so that if W is the energy 

21 

stored in the plasma oscillations before release, W ~ 10 ergs. 

2 

The number N of participating electrons is N = 4rrR 0 l 0 n, and if 

q 2Q 

R 0 = 2.5 x 10 cm, N ~ 6 x 10 . Thus each electron must possess, 

-9 

on the average, a energy w of 1.6 x 10 ergs. But if 6r cs is 
the charge separation associated with the oscillations and E is 



the electric field strength, we have from E = en5r , w = eE6r 

cs cs 

2 

E /n. Therefore E = 0.42 esu. 

We may now use Faraday's law to find the time rate of 
change of the magnetic field by assuming various values for the 
scale 6r E of the electric field. One possibility is to set 

~ lo > but f° r 3 large-scale magnetic field, possibly dipole- 
like, we might also set 5r ~ R . In any case, B ~ cE/er^, and 

Ij O Xj 

3*6 

for 6r = 6 x 10 cm, B ~ 2 x 10 gauss/sec, which implies a 
E 

rather large field if 6 ~ B, for a radio pulsar oscillation period 

8 • 

of about one second. However, for 6r E = 2 x 10 , B ~ 60 gauss/ 
sec, which is more reasonable. Note that if our estimate of 1 0 
is too low, then smaller values of B will suffice to produce 
the emission. 

It is observed that the pulses have characteristically short 
rise times. We note that mode coupling from the plasma 
oscillations to radiative plasma waves is non-linear and there- 
fore short groups of bursts with rather sharply defined frequency 
spectra might be expected to occur. 

Gyroradiation might also be a contributing source. For 
100 Me, the appropriate field strength is 36 gauss, but for non- 

relativistic electrons, the decay time for the electron kinetic 
17 

energy in such a field is roughly 5 x 10^ seconds. Thus, while 
it might be possible to get enough energy stored as electron 
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kinetic energy, many more orders of magnitude of participating 

electrons would be required for the necessary emissivity. 

However, gyroradiation can still contribute since fast rise- 

1R 

time bursts might result from stimulated "maser" emission 
operating in layers where the electron distribution is non- 
thermal. This would give nearly monochromatic emission, the 
band-width being determined by inhomogeneities in the magnetic 
field and by the electron relativity parameter. If the 
pulsation amplitude is large, the two radio bursts per stellar 
oscillation may have different characteristics. On the other 
hand, it is likely that the amplitudes are small since the 
calculations described in the previous section show that 
collapse occurs when the amplitudes are too large. 
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Table 1. WHITE DWARF MODELS WITH DEGENERATE ELECTRON PRESSURE ONLY 



cgs units 



Loq Density (qm/cc) 
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FIGURE CAPTIONS 


Outer radius (in units of 100 km) and central density 

6 

of a white dwarf (in units of 10 gm/cc) vs. time (in 

„ 32 

sec.). M = 8.38 x 10 g. 

Outer radius (in km.) of a neutron star as a function 


of time (in sec . ) . 



Pulsation Period = 1.335 sec 



t ( sec ) 
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